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Abstract 

Using the generalized Kolmogorov-Feller equation with long-range interaction, 
we obtain kinetic equations with fractional derivatives with respect to coordinates. 
The method of successive approximations with the averaging with respect to fast 
variable is used. The main assumption is that the correlator of probability densities 
of particles to make a step has a power-law dependence. As a result, we obtain 
Fokker-Planck equation with fractional coordinate derivative of order 1 < a < 2. 
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1 Introduction 



In studying of the processes with fractal time and long-term memory a generalized kinetic 
equation was proposed in [Ij. While the equation was of the master-type, its main property 
was the presence of the power-type kernel for a probability density to make a step. This 
type of equation was compared to the Kolmogorov- Feller equation in [2]. In this paper, 
we would like to go farther and to show the conditions under which one can obtain the 
fractional generalization of the Fokker-Planck equation from the Fokker-Planck equation. 

Fractional calculus [Sj HJ [5] has found many applications in recent studies in mechanics 
and physics, and the interest in fractional equations has been growing continually during 
the last years. Fractional Fokker-Planck equations with coordinate and time derivatives 
of non-integer order has been suggested in [6j. The solutions and properties of these 
equations are described in Refs. [2l |7j. The Fokker-Planck equation with fractional 
coordinate derivatives was also considered in [H [H HOl [HI [12] . 

The Kolmogorov-Feller equation is integro-differential one and it belongs to the type 
of master equations broadly used in different physical applications. It is well-known that 
the Kolmogorov-Feller equation can lead us to the Fokker-Planck equation under 
some conditions. In this paper, we use the method of successive approximations with the 
averaging with respect to the fast variable [Hj. The main assumption is that the correlator 
of probability densities of particle to make a step has power-law dependence. As a result, 
we obtain Fokker-Planck equations with fractional derivatives of order 1 < a < 2. 

In Sec. 2, the Kolmogorov-Feller equation for one-dimensional case is considered. In 
Sec. 3, we present a generalization of the Kolmogorov-Feller equation for two-dimensional 
case. The method of successive approximations is used for this generalized equation in 
Sec. 4. In Sec. 5, we use averaging with respect to the fast variable to derive fractional 
Fokker-Planck equations. Finally, a short conclusion is given in Sec. 6. 
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2 Kolmogorov- Feller equation for one-dimensional case 

2.1 Operator representation of the KF-equation 

Let P{t, x) be a probability density to find a particle at x at time instant t. The normal- 
ization condition for P{t,x) is 



/ + 00 
dxP{t,x) = l (t>0). 
-oo 



The Kolmogorov-Feller (KF) equation has the form 
dP{t,x) 



I dx'w{x')[P{t,x-x') - P{t,x)], P{0,x)=6{x), (1) 
where w{x') is probability density of particle to make a step of the length x', and 

+ 00 



dx'w{x') = l. (2) 

Let us introduce the operator representation of the KF-equation. We define the trans- 
lation operator 

T^' = exp{-x'dx}, (3) 

such that 

T,,P{t,x)=Pit,x-x'), (4) 
and the finite difference operator 

Ax' = /-T,', (5) 

where I is an identity operator. Then the Kolmogorov-Feller equation ([1]) can be presented 
as 

^^^ = L{A)P{t,x). (6) 
Here we use the integro-differential operator 

r+oo 

L(A) = - / dx'w{x')A^,. (7) 



The operator ([7]) will be called the Kolmogorov-Feller operator. 



2.2 KF-equation with fractional coordinate derivative 

The well-known fractional Caputo derivative of order a is defined by 

The fractional Marchaud derivative ^ of order a is defined by 

D:P(x) = ^|jP(.)-P(x)]^^^^, (0<a<l). (9) 
Using x' = X — equation ([9]) has the form 

If the function w{x') in KF-equation ([T]) is the exponential function 

MxO = ^^(^'), (10) 

where if(x') is a Heaviside step function, then Eq. ([1]) can be presented through the 
fractional coordinate derivative 

^^^ = aD^P(t,x), P(0,a;)=5(x), (0 < a < 1). (11) 

This is the Kolmogorov-Feller equation with fractional coordinate derivatives of order 
< a < 1. Note that the function w[x) is a probability density and it should sat- 
isfy the normalization condition ([2]). Then equation ffTOj) can be considered only as an 
approximation for large x' . 

2.3 Generalized KF-equation 

In the general case, we can suppose that probability density of particle to make a step of 
the length x' depends on the time instant t and the coordinate x. Then we should replace 
w(x') by w(t,x|x') in KF-equation ([T]). As a result, we can consider the equation 



di 



hoo 



dPit,x) ^. cia;'w;(t,x|a;')[P(t,a;-a;') -P(t,x)], P{0,x)=6{x), (12) 



oo 



where e is a small parameter. Here w{t,x\x') is the probability density to make a step of 
the length x' at the time instant t from the coordinate x. If w{t,x\x') = w{x'), then Eq. 
( fT2|) gives Eq. ([1]). We can assume that during any interval of time (t,t + dt) the value 
of the variable x{t) remains equal to x with probability 1 — p(t, x)dt and may undergo a 
change only with probability p{t,x)dt (see more about this equation in Sec. 55. of 
Then 

w(t,x\x') = p{t,x)w{x'), (13) 
where p{t, x) is a bounded function. If p{t, x) = 1, then equations (IT^ and ffT^ gives Eq. 

Using the operator ([5]), equation (fT2|) can be presented in the operator form 

dPit,x) 



dt 

where 



eL{t,x,A)P{t,x), (14) 



/ + 00 
(ia;'ti;(t,a;|a;') A^.. (15) 
'OO 

Equation (|T^ will be called a generalized KF-equation for one-dimensional case. 
2.4 Successive approximations 

The generalized Kolmogorov-Feller equation (IT^ can be rewritten in the integral form 

P{t,x) - P{0,x) = [ dr L{t,x,A)P{T,x). (16) 
Jo 

This equation can be presented as the integral Volterra type equation 

P{t,x) = P{0,x) +e [ dtiL{ti,x,A)P{ti,x). (17) 

Jo 

Let us consider the successive approximations. Substitution of equation (IT7|) in the form 

P{ti,x)=P{0,x)+e [ dt2L{t2,xA)P{t2,x) (18) 

Jo 
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into equation (|T71) gives 



P{t,x) = P{to,x) + £ [ dtiL{ti,x,A)P{0,x) + 



+ / dti / dt2L{ti,x,A)L{t2,x,A)P{t2,x). 



(19) 



Jo 



Changing the variables ti ^ t2, and t — > ti in equation (fTSll . and substituting into (fT9l) . 
we obtain an equation up to e"^ in the form 



If the function w{x') is the exponential function fllOp . then L{t,x,A) is a differential 
operator of order < a < 1 with respect to x, and L(ti,x, A)L(t2,x, A)P(t2,x) is a 
differential operator of the order < 2a < 2. To obtain fractional kinetic equations of the 
order < 2a < 2, we should consider the average procedure before a partial differentiation 
of Eq. fl2U]) with respect to time is realized. Without averaging, we derive equation of 
order < a < 1. 

3 Distribution function and Kolmogorov- Feller equa- 
tion for two-dimensional case 

Let P{t, X, y) be a function of probability distribution in a phase space. The variables x 
and y describe the phase space of a system. There are the following interpretations for 
the variables x and y. 

• A system with one degree of freedom can be presented by momentum x = p, and 
coordinates y = q. 

• A system can be described by action x = I, and phase y = 6. 



P{t,x) = P{0,x) + e / dtiL{ti,x,A)P{0,x)+ 





(20) 
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• ra-particle system can be defined by x = and y = (52, P2, (ln,Pn), or x = qi, 
and y = (pi, ^2,^2, gn,Pn). 

• The variables x describe a system, and y describes an environment of this system. 

We plan to use the reduced distribution and average values with respect to y, where 
y will be considered as a fast variable. 

3.1 Generalized KF-equation for two-dimensional case 

We assume that x e X cR and y eY cR, then r = {x, y) e X x Y C R'^ . We plan to 
consider x is a slow variable and y will be considered as fast variable. The distribution 
function P{t, r) in the region XxY satisfies the generalized Kolmogorov- Feller equation 



f-"- JXxY 

where d'^ri = dxidyi, and e is a small parameter. Here w{t, r|ri) is the probability density 
to make a step on the vector ri at the time instant t from the point r. 
Equation fl2T|) can be presented in the form 




(21) 



d_ 

di 



P{t,r) = eL{t,r, A)P{t,v), 



(22) 



where L{t, r, A) is a Kolmogorov- Feller operator 




(23) 



Here A^ is a finite difference operator 



A 



where Tr^ is a translation operator in X x y that is defined by 



= exp{-riV} = exp{-xi(9^. - yidy}. 
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We can assume that during any interval of time (t, t + dt) the value of the variable r(t) 
remains equal to r with probability 1 — p{t, r)dt and may undergo a change only with 
probability j9(t, r)(it (see Sec. 55. in [13]). Then 

w{t, r\r') = p{t,r)w{r'), (24) 

where p{t, r) is a bounded function. 

We can use the variables x, y instead of r. Then equation (|2T1) for the distribution 
function P{t, x, y) will be presented in the form 
dP{t,x,y) 



e dxidyiw{t,x,y\xi,yi)[P{t,x - xi,y -yi) - P{t,x,y)]. (25) 

JxxY 

This equation can be rewritten as 

d 

—Pit, X, y) = eL{t, X, y, A,, Ay)P{t, x, y), (26) 
where L{t, x, y, A^, Ay) is a Kolmogorov-Feller operator 

L{t,x,y,Ax,Ay)= / dxidyiw{t,x,y\x',y')[T^^Ty^ - I]. 

JxxY 

This is the generalized Kolmogorov-Feller equation in the operator form. 



(27) 



4 Method of successive approximations 

The generalized Kolmogorov-Feller equation 

^P(t,r) = £L(t,r,A)P(t,r) (28) 

can be presented as the integral Volterra equation 

P{t,T)=Po{r) + e [ dtiL{h,r,A)P{h,r), (29) 
Jo 

where Po{r) = P(0,r). Substitution of equation (129!) in the form 

P(ti, r) = Po(r) + £ / ' dt2 L{h, r, A)P(t2, r) (30) 







into Eq. (l29l) gives 

P{t,r)=Po{r)+e [ dh L{ti,r, A)Po{r) + 
Jo 

+ 6^ [ dti [ dt2L{ti,r,A)L{t2,r,A)P{t2,r). (31) 
Jo Jo 

Changing the variables ti ^ t2, and t — > ti in equation ( l30i) . and substituting into ( l3Ti) . 

we get 



'0 ^0 

Using the chronological multiplication 



P{t,r)=Poir)+e [ dh L{h,r, A)Poir) + 
Jo 

+ [ dh r dt2 L{h, r, A)L{t2, r, A)Po(r) + ... (32) 
Jo Jo 



TSn, AUr. AU J Hh,v,A)L{t2,r,A) h>t2; 
T{L{ti, r, A)L(t2, r, A)} = { (33) 

L(t2,r,A)L(ti,r,A) > ^i, 



equation fl32l) can be symmetric with respect to ti and t2: 

P{t,r)=Poir)+e [ dh L{ti,r, A)Po{r) + 
Jo 

+ l-e^ [ dh [ dt2T{L{h,r,A)L{t2,r,A)}Po{r) + ... (34) 
^ Jo Jo 

This is the symmetric representation of equation (l32l) . 



5 Averaging with respect to the fast variable 

In this section, let us consider the variables r = {x, y) as slow (x) and fast {y). Substitution 
of ^ into ([31D gives 

P{t,r)=Po{r)-6 [ dh [ d\iw{ti,r\ri)ArMr) + 

Jo JxxY 
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t ft 



+ -s' dh dt2 / rfVi / dV2T{w;(ti,r|ri)Ar,w(t2,r|r2)ArjPo(r) + 0(£^). 

^ Jo Jo JxxY JxxY 

(35) 



The function w{t, r\ri) is the probabihty density to make a step of the vector ri at 
the time instant t from the point r. The first assumption is that this probabihty density 
has a weak dependence (up to terms of order 0(e)) on the point r, i.e., 

w{t, r + nlra) = w{t, rjra) + 0{e). (36) 

Then 

An w{t2,r\v2) = 0, 

and 

w{ti, r|ri) An w{t2, r|r2) A^ = w{ti, r|ri) w{t2, r|r2) A^ A^^. 
As a result, equation (1551) has the form 

P{t,r)=Po{r)-e [ dti [ d\,w{h,r\ri) A,Mr) + 

Jo JxxY 



+ le^ [ dti [ dt2 [ d\i [ d\2T{w{ti,r\ri)w{t2,r\r2)}Ar,Ar,Po{r) + 0{e^). 

^ Jo Jo JxxY JxxY 

The second assumption states that Po{r) = Pq^x^u) does not depend on the fast 
variable y up to e-terms such that 

Po{r) = P,{x,y)=p^{x) + 0{e). (37) 

Then, we have 

An Ar, Po(r) = A,, A,, po{x) +0{e). (38) 
Averaging over the variable y will be denoted by < >y. We also use the notations 

p(t,x) =<P(t,r) >„ (39) 

and p(0, x) = po{x). 

Averaging of equation (!35|) with respect to the fast variable we obtain 

p{t,x) = po{x) - e dti / dxiA{ti,x\xi) Ax-^po{x) + 
Jo J X 
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+ / dti I dxi I dx2B(ti,x\xi,X2) A^^A^^po{x) + 0{e 
^ Jo J X J X 



3\ 



(40) 



where we introduce the functions 



A{ti,x\xi) = j dyi < w{ti,r\ri) >y, (41) 

B{ti,x\xi,X2) = / dt2 / dyi / dy2 < T{w{ti,r\ri) w{t2,r\r2)} >y . (42) 
Jo Jy Jy 

Using r = {x,y), these functions can be presented by 

A{ti,x\xi) = j dyi < w{ti,x,y\xi,yi) >y, (43) 

B{ti,x\xi,X2) = / dt2 / dyi / dy2<T{w{ti,x,y\xi,yi)w{t2,x,y\x2,y2)}>y- (44) 
Jo Jy Jy 

The third assumption is that the function B{ti,x\xi,X2) is diagonal with respect to 
variables xi and X2 up to e-term, i.e., 

B(ti,x\xi,X2) = B(ti,x\xi)6{xi — X2) + 0{e). (45) 

Then the operator A^.^ A^,^ will be the finite difference operator A^^ of second order. 
This allows us to have fractional derivative for the exponential function B{t,x\xi) since 
Marchaud and Riesz fractional derivatives [3] are defined through the finite difference 
operator. 

The fourth assumption is that the functions ^4(^1, a;|xi) and B{ti, x\xi) are exponential 
functions up to e such that 

A{t, x\xi) = ait) / 1+1 -^(^1) + 0{e), (0 < ai < 1), (46) 

B{t, x|xi) = b{t, x) / \ h{xi) + 0{e), (1< as < 2), (47) 

K[a2, ^) \Xi\ ^ 

where H{x) is the Heaviside step function, and 



«:(«, n) = -r(a)A4a), A^{a) = E^"^)''' fc!(n - fc)! 
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with n — 1 < a < n. 

As a result, Eq. fj40l) gives 

p{t, x) = po{x)+e [ dti a(ti)D°Vo(a;) + [ dh K^u x)B'^^po{x) + 0{e^). (49) 
Jo Jo 

Here is Marchaud fractional derivative [3j of order a with respect to x, which is defined 
by the equation 

1 r<^ A" f(x) 

°°^w = ;<^i <"<«)■ (50) 

where is a finite difference of order n such that 

n . 
^ r 71 ' 

A^/(^) = {1- TyTfix) = ^(-ir „ • - "^2/)- (51) 

— ^ m! n — m ! 

m=0 ^ ' 

In general, the variables x and Xi can be vectors in the A^-dimensional space R^, 
where = 1, 2, 3, .... The fourth assumption for the functions Ait^, x\x\) and -B(ti, x\x\) 
can be realized in the form other than ( H6l) and ( H7I) . We can suppose that the functions 
A{t\^x\x\) and -B(ti,x|xi) are exponential functions up to e such that 

A(t,x|xi) = a(t);p7^ v i iLa, + ^(^)' (0 < «i < 1), (52) 



x|xi) = a:) , .^^^^ + 0(e), (1< as < 2), (53) 

dAr(2,a2) |a;i|^+"2 



where x, Xi G M^, and 



2-"7r^+^/M„.(a) 
sin(«7r/2)r(l + «/2)r((Ar + «)/2) ^ ^ 



with n — \ < a < n. As a result, Eq. fHUl) gives fractional equation of the form 
where D° is the fractional Riesz derivative (see [3J Sec. 25. 4) of order a with respect to 
X G R^, which is defined by the equation 

D°/(x) = , / , / d^'y, {n-\<a<n). (55) 
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To denote Riesz fractional derivative ( |55i) is also used the notation (— A)°/^. Note that 
the Fourier transform T of this derivative (see Property 2.34 of [4J) is defined by 

(-^{D:/(^)})(k) = |kr (w(x)})(k). 

The representation of the assumption in the form fl52l) and fl53l) allows us to obtain 
fractional kinetic equation for arbitrary A^-dimensional space (for example, in the 3- 
dimensional Euclidean space). 

The partial time differentiation of equation (jl9|) gives 

^p(t, X) = £a(t)D^Vo(a;) + e^b{t, x)DTpo{x) + 0{e^). (56) 

Substitution of equation (149!) in the form 

Po{x)=p{t,x)-6 [ dha{h)B'^'po{x) + 0{e^) (57) 
Jo 

into equation (156!) gives 

—pit, x) = sa{t)-D'^' p{t, x) + (b{t, x)D^2 - c(t)Df i)p(t, x) + 0(5^), (58) 

where 

c(t) = a{t) / dti a{ti). 
Jo 

Equation (158|) up to 0{6^) has the form 

=e«WD;V(*..) +e=(KM)Dr -e(*)D-),(*..). (59) 
This is the fractional kinetic equation with noninteger derivatives of the order 1 < 0^2 < 2 
and < 2ai < 2 with respect to coordinate x. 

If a{t) = 0, i.e.. A{t, x\xi) = 0, then we have the fractional equation of order 1 < 02 < 
2 with respect to x such that 

^^p{t,x)=e'b{t,x)T>Tp{t,x). (60) 

This is the fractional Fokker-Planck equation that is suggested in [6] to describe fractional 
kinetics. The solutions and properties of these equations are described in Refs. [21 [7]. 
The Fokker-Planck equation with fractional coordinate derivatives are also considered in 
l8l[9l[Tni[IIl[T2]. 
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6 Conclusion 



In this paper, the Fokker-Planck equations with coordinate derivatives of non-integer order 
1 < a < 2 are derived by using Kolmogorov-Feller equation. The method of successive 
approximations and the averaging with respect to a fast variable are used. The main 
assumption in the paper is that the correlator of two probability densities of particle 
to make a step has the form of exponential function. As a result, the Fokker-Planck 
equations can have derivatives of order 1 < a < 2. 



References 

[1] E.W. Montroll, M.F. Shlcsinger, "On the wonderful world of random walks" in: 

J. Lebowitz and E.W. Montroll, Editors, Studies in Statistical Mechanics, Vol. 11, 

North-Holland, Amsterdam (1984), p.l. 
[2] A.I. Saichev, G.M. Zaslavsky, "Fractional kinetic equations: solutions and apphca- 

tions" Chaos 7(4) (1997) 753-764. 
[3] S.G. Samko, A. A. Kilbas, O.I. Marichev, Fractional Integrals and Derivatives Theory 

and Applications (Gordon and Breach, New York, 1993). 
[4] A. A. Kilbas, H.M. Srivastava, J.J. Trujillo, Theory and Application of Fractional 

Differential Equations (Elsevier, Amsterdam, 2006) 
[5] I. Podlubny, Fractional Differential Equations, (Academic Press, San Diego, 1999). 
[6] G.M. Zaslavsky, "Fractional kinetic equation for Hamiltonian chaos" Physica D 76 

(1994) 110-122. 

[7] G.M. Zaslavsky, "Chaos, fractional kinetics, and anomalous transport" Phys. Rep. 
371 (2002) 461-580. 

[8] G.M. Zaslavsky, " Renormalization group theory of anomalous transport in systems 
with Hamiltonian chaos" Chaos 4 (1994) 25-33. 



14 



[9] A.V. Milovanov, "Stochastic dynamics from the fractional Fokker-Planck- 

Kolmogorov equation: Large-scale behavior of the turbulent transport coefficient" 

Phys. Rev. E 63 (2001) 047301. 
[10] V.V. Yanovsky, A.V. Chechkin, D. Schertzer, A.V. Tur, "Levy anomalous diffusion 

and fractional Fokker-Planck equation" Physica A 282 (2000) 13-34. 
[11] V.E. Tarasov, "Fractional Fokker-Planck equation for fractal media" Chaos 15 (2005) 

023102. 

[12] V.N. Kolokoltsov, "Generalized continuous-time random walks (CTRW), subordina- 
tion by hitting times and fractional dynamics" E-print larXiv:0706. 19281 

[13] B.V. Gnedenko, The Theory of Probability (Chelsea, New York, 1962). 

[14] G.M. Zaslavsky, Chaos in Dynamic Systems, (Harwood Academic Publishers, New 
York, 1985) 



15 



